THE HARTMAN- WATSON DISTRIBUTION REVISITED: 
ASYMPTOTICS FOR PRICING ASIAN OPTIONS 



STEFAN GERHOLD 



Abstract. Barrieu, Rouault, and Yor [J. Appl. Probab. 41 (2004)] determined 
asymptotics for the logarithm of the distribution function of the Hartman- 
Watson distribution. We determine the asymptotics of the density. This re- 
finement can be applied to the pricing of Asian options in the Black-Scholcs 
model. 



1. Introduction and main result 

The distribution of the integral of geometric Brownian motion has attracted a 
lot of interest, in particular because it is needed to calculate the price of Asian 
options in the Black-Scholes model. Yor [16] found the formula 

(1) P[A^ e du \ W t + vt = x] = — exp - - — I (e x /u)f e * /u (t) du, 

u \ It 2u J ' 

where I„ denotes, as usual, the modified Bessel function of the first kind, and 

ri 



A| v) = f exp(2(W h +fh))dh, 
Jo 



where W is a standard Brownian motion. The present note focuses on the func- 
tion f r (t) in ([I]) , which is the density of the Hartman- Watson distribution [15] . 
It is defined for a positive parameter r by the Laplace transform 

e- ut f r (t)dt= 1 ^±, R( u )>0. 
o W) 

Small time asymptotics of the conditional density ([T]) correspond to left tail asymp- 
totics of f r (t). Numerical problems in the evaluation of ([TJ for small t prompted 
Barrieu, Rouault, and Yor [1] to analyze the left tail of the Hartman- Watson dis- 
tribution asymptotically. Using the Gartner-Ellis theorem from large deviations 
theory, they obtained the asymptotics 

(2) F r (t)=exp( ^^+o( S Y '- )j, t^O, 

for the distribution function. However, this result is not immediately applicable 
to the calculation of ([T]). Barrieu et al. [I] write that "the standard asymptotic 
methods (e.g. the saddle point method) do not seem to be suitable for this study" , 
and that they "are not able to refine these results for the Hartman- Watson density 
itself". 



Date: May 9, 2011. 

2010 Mathematics Subject Classification. Primary: 62E20; Secondary: 60J65. 

Key words and phrases. Hartman- Watson distribution, Asian option, saddle point method. 



1 



2 



STEFAN GERHOLD 



In fact the saddle point method [14] is applicable to the Laplace inversion 
integral 



+ioo T — ( r \ 



e 



du 



Mr) 

for the density f r (t), but needs some care. First, replacing Iy^j(r) by an asymptotic 
approximation relieves us from studying potentially difficult monotonicity proper- 
ties of the modified Bessel function, and allows to formulate the result in a way 
that avoids roots of equations involving the Bessel function. Second, it turns out 
that elementary approximations of the integrand's saddle point lead to integration 
contours that are too far away from the saddle to make the method work. An ap- 
proach based on a contour through the exact saddle point establishes the following 
asymptotics for the density f r (t). For brevity, we write 

p = log — 

Theorem 1. For t > 0, denote by u Q (t) the largest solution of the equation 

logu p 1 

3 t = —7=- ~t= + ~r> 

which exists for all sufficiently small t. Then the Hartman- Watson density satisfies 



(4) f A t) = ^ log ^ 2 _ 2p x (l + 0(VHog(lA) 2) ) 

(5) _ y/e log(lA) 1 / 2 tMt)+ ^r } ( / loglog(lA) \ 

(5) -27TlnM t C 1 + V Wl/t) ) 



as t — > 0. 



Formula Q gives a much better approximation than ([5]) ; the simplification in ([5]) 
is of little use, since uo{t) has to be computed anyway to evaluate ((U) or (J5J) nu- 
merically. 

To get a feel for the growth of the exponential in Q, we expand uo(t) by boot- 
strapping (cf. de Bruijn [2, Section 2.4]): 

ffil „ ^- 1Q S(1A) 2 (-, 1 21oglog(l/t) 2p + log2 / 1 
Therefore the exponent in Q has the expansion 

(7) -^^.-M^-S^W) 



2t t 

t V t 



(l + P + ilog2)^ + of^). 



This shows in particular that the formula 

( log(iA) 2 /log(iA) 2 

(8) / r (t)=exp(^ ^L^+o[ B y ; 

obtained from ([5]) by formal differentiation, is correct. 

For numerical accuracy, it is certainly preferable to use Q as it is, without 
replacing the exponent by ([TJ); still, the expansion (j6|) can serve as good initial 
guess when computing the root of ((3|). In this way, the leading term of f r (t) can 
be calculated effortlessly even for extremely small values of t, say t = 10~ 50 . 
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2. Analysis of the Laplace inversion integral 

The Laplace inversion formula yields the representation 

i nR+ioo t — ( r \ 

fr{t) = rr- / e du, 
2itt Jij-ioo Io(r) 

where i? > 0, so that the integration contour lies in the right half-plane. To 
estimate the growth of f r (t) near t = 0, we have to investigate the singularity at 
infinity of the integrand. For large index, the modified Bessel function admits the 
expansion [51 [T3] 

(9) L(r)~(0V^- 1 /2( Co + 2l + ...), 

where the Cj are constants, with cq = l/v2~7r. This holds for v — > oo, uniformly 
w.r.t. arg(^), as long as arg(f) is bounded away from ±7r. Horn j8j shows (O 
for J, y (r), the Bessel function of the first kind, but from the relation l v (r) — 
e _1 ""/ 2 J ! ,(re 7rl / 2 ) one easily sees that replacing J„(r) by L(r) only affects the con- 
stants c\ , C2, . . . in this expansion. If we let the real part R of our integration 
contour tend to infinity as t — ¥ 0, we therefore have 

f*® = to W2-T < \ P +1Q ° (^^^(V^r^-^du x (i + ocjt 1 ")) 

(27r)- i / 2 il (r) J fl _ ioo V27 

2-7/4 pR+ioo 

= „ /0 . T , s / exp (ttt - i v / 2lxlogu + V2(l + 
7T J / 2 iI (r) 7 fl ,_ ioo V 

- Jlogu) dux (l + CKiT 1 / 2 )). 

The integrand of the latter integral has a saddle point, let us call it uq = uo(f ), which 
is found by equating the derivative to zero. This yields equation ([3]). Shifting the 
integration contour through the saddle point achieves concentration, so that only 
a small part of the contour matters asymptotically. 

In many instances of the saddle point method, it suffices to choose a contour that 
passes through an approximation of the saddle point. In our example, one might 
try to use a contour based on the first terms of ©• However, painful calculations 
reveal that the concentration of the integrand around the approximate saddle point 
is insufficient, no matter how many terms of ([6]) are taken. We therefore set the 
real part of the integration contour to the exact saddle point, so that R = uq: 

2-7/4 /-uo+ioo _ . 

(10) f r (t) ~ 3 / exp(^-iv / 2Mlogw+v / 2(l+p)v / ^-|logw)dw. 

Let y denote the new (real) integration variable: 

u = uq + iy, — oo < y < oo. 

Close to the saddle point, i.e. for small values of the new integration variable y, we 
have the uniform expansions 
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and 



r-i , — , 12/ . ilog(u )y log(u )y 2 /log(u )y 3 



We insert these into the exponent of (fTOj) and obtain 

ut — ^V^wlogw + \/2(l + p)\fu — -j logw 

= u t- i\/2u^logMo + \/2(l + p)\M)- iloguo 
'log(u )y 3 ^ 



(11) -Afy 2 + 0(- 



5/2 



where 

(13 ) - t ( 1 1 pf jggjgjft/g 

1 J -21og(lA)2l v 1 + U l log(lA) 

Note that the y-terms in (jTTJ) vanish, because we integrate through a saddle point. 
We now have to identify a range 

—h <y < h 

for y — 3(u) that captures the main contribution to the integral (fTQ|) . A good 
choice is 

■ _ log(iA) 2 

t S/2 > 

because it satisfies h\M — > 00, so that the integral of the local expansion (fTTj) can 
be completed to a full Gaussian integral: 



h 1 r hV2M 

e" 

h \/2M J-hV2M 



1 r hV2M 



1 /" 

(14) - -= / e-^^dw 

tT V^rlogCl/*) 



M i 3 / 2 

Moreover, the error from (fTTj) . the local expansion at the saddle point, is o(l), since 



(15) ii^=0 ( Vtlog(lA) 2 ). 

We can thus determine the asymptotics of the portion \$s(u)\ < h of the inte- 
gral ([10]): 

2-7/4 ruo+ih _ 

( 16 ) o**/o * t ; n / exp ( iti — i y/2u log it + \/2( 1 + p)y/u — \ log u ) dtt 
Tr d / 2 il (r) y„ _ ih V J 

q — 7/4 1 1 r-h 

a u t- 2 V2"o iog« +v'2(i+p)\/uo- 4 iog«o / e~ My2 (iy 



TT 3 / 2 1 (r) 



-A 



o-7/4 1 
/i7^ 71 f— 1/2. — 1/4 „«ot— «V5 log ti + V^(l+p) v /Mo 

(17j 
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This gives the right-hand side of (jlj, after expressing ^/uologuo via the saddle 
point equation ([3]), which yields 

(18) - l\/2u logw = ~2u t - py / 2u + |. 

Furthermore, expanding uq by ([6]) gives the expression in (J5J). Note that we have 
not yet proved ((4]) and ((5]); it remains to show that the tails of (flO|). i.e. the parts 
where |3(u)| > h, are asymptotically negligible. This ensures that (TIB]) indeed 
captures the asymptotics of / r (t). 

3. Tail estimate 

To bound the tails of (flU)) . it suffices to consider the case y = > h, since 

the lower half of the tail can be handled by symmetry. We first deal with the part 
of the contour in (|10j) where the imaginary part of the integration variable is very 
large, say y > e log< - 1 /*- 1 / 4 . Then y clearly dominates uq, and it follows from 

K(logu) ~ logy, 3?(\/u) ~ and Ji^logw) - ±A/2ylog|/ 

that the absolute value of the integrand is bounded by 

for small t. Hence we obtain the bound 

(19) e" ot / e-^d?/~2exp( M oi+slog(lA) 2 -e los(1/t)2/8 ) . 

J G lo g (l/t)2/4 V / 

Finally, we bound the portion of the integral flT0|) that is close to the central 
part, i.e., 

(20) h < y < fpxW*)'/*. 

The following lemma shows that, for small t, the absolute value of the integrand 
decreases as y increases. 

Lemma 2. Let B be a real number. Then, for $t(u) > and \u\ sufficiently large, 
the real part of y/ulogu + B ^/u decreases w.r.t. \$s(u)\. 

Proof. We write u — x + ly. By symmetry, it suffices to consider the case y > 0, so 
that arg(it) > 0. Straightforward calculations show that 

j^(v^) = |M 3/2 (ycos^^i -xsinHMHi) 

and 

|-5R(v^log u) =\ M 3 / 2 ^ (log M + 2) (y cos ^ - x s i n 2^*1) 

— arg(it) I x cos ^ v ' + y sin ^ J 1 I. 
Hence we are led to investigate the sign of 

y ((log M + B + 2) (cos ^ _ 2 sin «s£l) _ arg ( n ) (2 CO s ^ + sin hMh! ) ) . 
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Suppose that |u| is so large that log \u\ + B + 2 > 12. In the preceding formula, we 
estimate the trigonometric functions by the first term of their Taylor series at zero, 
except the first cos, where we use two terms. This yields the lower bound 

V (l2 (l - \ arg(u) 2 - ^ arg(u)) - | arg(u) - \ arg(w) 2 ) 

= y(l2-2arg( W ) 2 -^arg( W ) 

= y (12 — 2 arctan(w) 2 — ^ arctan(u;)) \ w _ y , x ■ 

Now observe that arctan(iu) 2 < 7r 2 /4 and arctan(w) < w for w > 0, so that 

12-2 arctan(w) 2 - ^ arctan(w) > 12 - tt 2 /2 - 7 > 0. 

This shows that 9fi(y / ulogu + B^/u) has a positive derivative w.r.t. y. □ 

Therefore, we can bound the part (|20|) of the integral (flOl) by the value of the 
integrand at y — h times the length of the path. By (fTTj) . ([T8|l . and 

M y 2 Ufc~i lo g(i/*) 2 , 

this amounts to a bound of the form 



(21) e -*«o(t) + %/2(l+2p)VMt)-^log(l/t) 2 +o(log(l/t) 2 ) xe ilog(l/t) 2 



_ e -t«o(t) + V^(l+2p)^« (t)-ilog(l/t) 2 +o(log(l/t) 2 )^ 

To complete the proof of Theorem [1] let us now compare the six error terms that 
arose in the analysis (see Table [1]). Note that the error from completing the tails 
of the Gaussian integral in (fl4)l is 



lhV2M V M W 

(22) = e~i 

If M is not expanded, i.e., (fT2|) is used, then the error from the local expansion 



w=hV2M 
±log(l/t) 2 +o(log(l/i) 2 ) 



Source of error 


Relative error 


Replace I v by ((9|) 


0(f/log(l/f)) 


Local expansion (see (HiT) and (fT5l)) 


O(^log(l/t) 2 ) 


Gaussian tails (see (fl4l) and (l22l) 


exp(-ilog(lA) 2 +o(log(l/t) 2 )) 


Relative error of M (see ((131) and (|14|) 


n /log log(l/t) \ 
U( - l°g(l/t) ' 


Outer tail (see ([19])) 


CX P (-C 1 °g( 1 /t) 2 /8 + ( e log(l/t) 2 /8)) 


Inner tail (see ([2T])) 


exp(-ilog(lA) 2 +o(log(lA) 2 )) 



Table 1. 



dominates, which leads to (U]). If, on the other hand, the expansion (| 13[) of M is 
taken, then it is the relative error of M that prevails. 
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4. Comments 

The left tail of the Hartman- Watson distribution (see ([8])) is somewhat thinner 
than that of the Levy distribution (stable distribution with index a = i), with 
density 

g(t) = ^=e- x ' 2 \ t>0, 

and Laplace transform 

poo 

/ e~ ut g{t)dt = e-^, > 0. 

Jo 

The faster decay of the Laplace transform of the Hartman- Watson distribution, of 
order exp(— y/u log u), becomes manifest in the additional factor log(l/i) 2 in the 
exponent of flSJ). 

We now briefly comment on possible refinements of Theorem [TJ Technically 
speaking, continuing the expansion (fi"3]) and inserting into ([IT]) refines J5} to a 
full asymptotic expansion. A better expansion, respecting the asymptotic scale of 
the problem, can be obtained by retaining the explicit formula (|12[) for M, and 
taking more terms in (j9]) and (fTT]) . This should pose no essential difficulties; note, 
however, that each term in the expansion Q gives rise to a new saddle point, as the 
coefficient of 1/u in (J3j> changes. Thus the expansion will involve several implicitly 
defined functions of t besides u (t). Since the dependence of the solution of © on 
the coefficient of 1/u is light, it might be possible to give an expansion that features 
only Uo(t). This seems of little practical interest, though. 

Concerning applications, our results can be used as a substitute for the Hartman- 
Watson density for small arguments, in particular, for evaluating the density of A^ 
numerically for small t (after integrating (fTJ) w.r.t. the law of W t + vt). The related 

problem of determining small time asymptotics for the density of is left to 
future research. This density is difficult to evaluate numerically for small time [5]. 
Analyzing it asymptotically requires handling a double integral; see Tolmatz |ll[|12j 
for asymptotic evaluations of double integrals pertaining to other functionals of 
Brownian motion (resp. the Brownian bridge). For v = 0, the analysis should be 
simpler, as the density of A\ can be expressed as a single integral, via Bougerol's 
identity [H [10] . Note that tail asymptotics [6] and large time asymptotics (4j [10] 
of A[^ arc known. See also Dufresne [5] for related limit laws. 
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